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^) ' We present a general method for studying coupled qubits driven by adiabatically changing external 

parameters. Extended calculations are provided for a two-bit Hamiltonian whose eigenstates can 

be used as logical states for a quantum CNOT gate. From a numerical analysis of the stationary 

Schroedinger equation we find a set of parameters suitable for representing CNOT, while from 

a time-dependent study the conditions for adiabatic evolution are determined. Specializing to a 

concrete physical system involving SQUIDs, we determine reasonable parameters for experimental 

purposes. The dissipation for SQUIDs is discussed by fitting experimental data. The low dissipation 

obtained supports the idea that adiabatic operations could be performed on a time scale shorter 

2*t: , than the decoherence time. 
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I. INTRODUCTION 



The basic elements for processing quantum information (e.g. to perform quantum computation pj) are quantum 
bits (qubits), namely two level systems exhibiting quantum coherence between the states, quantum register (arrays of 
qubits) and quantum gates Q ■ Computations are performed by the creation of quantum superpositions of the qubits 
and by controlled entanglement of the information on the qubits . The main goal of any physical implementation of 
a quantum information-processing device is therefore to control systems of coupled qubit with a phase coherence time 
long enough to permit the necessary manipulations. Various physical systems have been proposed for the physical 
implementation of qubits .including photons [3J, trapped ions [aL spins in nuclear magnetic resonance (NMR) |(|, 
. £h ' electrodynamics cavities @ and semiconductor quantum dots [Sj, to mention some of the most popular. Other 
experiments and proposals focus on superconducting Josephson devices, where almost macroscopic (mesoscopic) 
devices like Josephson junctions [j|, SQUIDs or Cooper pair boxes [ll|. that is devices fabricated in condensed 
matter physics, are brought to behave quantum mechanically. These devices promise certain advantages like large- 
scale integration and fabrication as well as ease of integration with conventional electronics. To manipulate qubits 
quantum gates are necessary, that is logic devices capable of operating on linear combinations of input states. 

Among the possible mechanisms for manipulating coupled qubits, adiabatic procedures |12| are of special interest. 
Quantum adiabatic evolution provides a natural framework for solving combinatorial search problems on quantum 
computers [13] . Any problem which can be recast as the minimization of an energy function (which can then be 
converted into a quantum Hamiltonian) can potentially be solved by adiabatic quantum computation. General 
problems have been treated numerically, and studies of a set of Exact Cover instances designed to be hard have shown 
polynomial behaviour out to instances containing as many as twenty bits |l4j . Whereas a conventional quantum 
algorithm is implemented as a sequence of discrete unitary transformations that form a quantum circuit involving 
many energy levels of the computer, the adiabatic algorithm works by keeping the state of the quantum computer 
close to the instantaneous ground state of a Hamiltonian that varies continuously in time. Therefore, an imperfect 
quantum computer implementing a conventional quantum algorithm might experience different sorts of errors than an 
imperfect adiabatic quantum computer. In fact, an adiabatic quantum computer has an inherent robustness against 
errors that might enhance the usefulness of the adiabatic approach ^^|. Local operations on single qubits (such as 
NOT) or two coupled qubits ( such as the adiabatically controlled CNOT gate we shall discuss) [17|][lg] are a ^ so 
possible, where adiabatic operations take place as a sequence of discrete transformations acting on a few qubits at 
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a time. In this line it is important to study a possible trade-off between the advantages of error reduction due to 
adiabatic evolution and the longer times required for gate operations. 

In this paper we will explain some general principles for studying adiabatic SQUID qubit operations focusing 
particluarly on a CNOT gate, and present numerical calculations relevant to the behavior and design of such systems. 



II. COUPLED-QUBIT HAMILTONIAN 



In discussing coherence properties of the SQUID under adiabatic inversion, we have suggested its interest for the 
elements of the "quantum computer" [l6l |. The single-bit NOT operation can be realized by adiabatic inversion. 
The next most complicated operation is the two-bit operation CNOT, with which a computer may, in principle, be 
constructed. CNOT is a two-qubit operation and we will try to represent it by two interacting double-potential well 
systems. Qualitatively, we will use the procedure of performing an adiabatic NOT on the first qubit while trying to 
influence its behavior by the state of the second. We'll find a region of parameter space where this works. 

In the one dimensional or one SQUID problem one has a Schroedinger equation in the variable <E> with a kinetic 
term and the potential term 

U= { - —J -^cos(27r*/$o). (1) 

which for small <!> yields a double well potential. I c and L are the Josephson critical current and the inductance of 
the superconducting ring respectively, while & ext is an applied external flux. When Q ext is swept slowly as a function 
of time, as explained in ref [la ], an adiabatic inversion or level crossing can be induced, amounting to a realization of 
the NOT operation. If the state is originally in the left potential well it is transferred to the right well and vice-versa, 
implying a reversal of the current direction in the superconducting ring. 

We now wish to investigate this idea of quantum gates generated by adiabatic transformations to systems of more 
than one variable, in particular for the two variable CNOT operation. Although we are only concerned here with two 
SQUIDS, we briefly indicate a method valid for many qubits. The equation for an array of underdamped flux-linked 
rf SQUIDs is 

$ - $ ext = Li (2) 

In this equation <E> and i are meant as column vectors representing all the fluxes $j and currents ij in the j-th 
rf SQUID loop, & ext is the column vector corresponding to the external fluxes, while L is a matrix representing 
the self- and mutual- inductances. The current ij in j-th ring is expressed in terms of the capacitance Cj and the 
superconducting Josephson current Ij : 

where $ = f§ = 2 • 1CU 7 G cm 2 is the superconducting flux quantum and R is the effective resistance of the junction. 

An important property of L, by the reciprocity of mutual inductances, is that L is a symmetric matrix. As we need 
i for the linear homogeneous system (LHS) of Eq[3j, we invert L: 

i = L- l {$-$ ext ) (4) 
Since L is symmetric, L^ 1 is also symmetric. Neglecting the dissipative term, we now have Eq [3] as 

- E L 7k(* ^ - = C & ( 5 ) 



k 



that can be written as 



by introducing the potential 



dU 



= Cj*i (6) 
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j-k 

Finally switching to the reduced dimensionless flux variable (j> = 3>fp the last two equations become 



= (f^ E L 7^ - 0*0* - - (f^E^coflfc (9) 

To make this situation look more symmetric, we follow recent practice and introduce the shifts <j> — > + tt, 
ea:t — > </> e:Et + 7r, which move the maximum of the coscj) term to — [l^. This does not affect the quadratic term, 
but it must be kept in mind that 4> ext = now corresponds to a non-zero applied field. We thus have finally 

U->U= {^?\ J2 Ljfa ~ <t> ext )k{<j> - cb ext ) 3 + (^Erfcoah , (10) 
where henceforth we use the shifted variables. 



A. Two variables 



We now specialize to two SQUIDS, called 1 and 2, with variables <pi, <f>2- In the absence of mutual inductance, each 
one can be thought of as a qubit, whose dynamics is described by a double- well potential. For CNOT we shall think 
of 4>\ as the target bit and <f>2 as the control bit. We shall apply a sweeping flux (j)\ xt and would like that this sweeping 
flux induce an adiabatic inversion if the control bit is |1) (e.g. current flowing clockwise in SQUID 2), and not induce 
this inversion if the control bit is |0) (current flowing in the opposite direction in SQUID 2). 

L\ L\2 

L\2 L2 



For the two variable system the matrix L 



can be inverted, giving 



L~ = 



1 



L1L2 



T 2 
^12 



Li 
—L12 



-L12 
Li 



(11) 



L12, the mutual inductance between the two rings, could be controlled experimentally by a further device with 
Josephson junctions |20| and switched off for convenience in performing other operations. 
Writing out U we obtain 



U = (^) 2 o77 -^[L 2 (0i-0rT+£iW2-^) 2 -2L 12 (^^^^ (12) 

Let us first look at this potential with the <p ext zero, and imagine varying <pi at fixed 4>2- For L12 — 0, we would 
simply have the usual symmetric double potential well for <f>\. Now as h\% is turned on, a tilt is introduced into this 
4>i potential. A coupling term ~ ^i2</ , i02 is added providing a bias in the potential for (f>i, so that the double well is 
asymmetric even though <pf = 0. The direction of this bias depends on whether $2 is positive or negative. 

If we think of $2 as essentially fixed, and now suppose sweeping (j>f xt , we see that this sweep will either increase 
the asymmetry present in the double well, or decrease it. If the wells are caused to be further separated, we have no 
inversion, if the wells are brought together and cross, we will have an inversion and a "NOT" . Which case occurs will 
depend on the sign of <f>2- The condition for CNOT is accomplished: according to the state of $2 , an inversion takes 
place or does not take place in the 4>x variable. 

The discussion is more convenient if we introduce dimensionless parameters 

1 VX1X2 2tt 2tt 

I L lL2 L\ 2 Pl ~W Lh fo-lF Ll2 (13) 
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as well as the dimensionless inductances 



h = -^L= h = ~^= h2 = ^= (14) 

\JL\Lii \J L1L2 v L\Lii 



We may then write Eq^21 as 



U = {^?\{\[H<t>i ~ tt Xt ) 2 + Ufa - 0D 2 - 2^12(02 - 0D(0i - + Picos^i + focosfo} (15) 

With this potential the Schroedinger equation for the wavefunction tjj((j)i,4>2,t) is 

iip = Hip (16) 

with the following Hamiltonian 

-1 d 2 -Id 2 

7-1 = T TT79+ J. TTJo+U 17) 

2d(&) 2 d<ft 2C 2 (^) 2 ^| 

Introducing 



C = VdC 2 (18) 
as the typical capacitance, and defining the energy scale 

E a = 1/VIC (19) 
H. can be cast in the form of an energy times a dimensionless Hamiltonian 

H = EqH . (20) 



H=—— + 7T- 7 ^ + V (21) 



-1 a 2 -\__d 2 

with V 

V = Voi^hfa - <P\ xt f + h(fo ~ (t>2 Xt f - 2^12(^2 - - + Picosfa + focosfo} , (22) 



and the dimensionless parameters 



Hi = C 1 E A 2 fl2 = c 2 E (% 2 Vo = ^CJlA 2 . (23) 

Ztt Ait Ait 

Fig. shows the equipotential contours of V((f>i, </> 2 ), with its four potential wells. 

We use natural units: h = l,c = l,e 2 /hc = 1/137, thus $ = hc/2e = 7r v / 137, and (§£) 2 = 137/4 = 34.3; also 
yj farad/ henry = 30. Using these values 



Note that these parameters are not all independent, as the following relation exists 

VJ 1 ^ = Vo ■ (25) 

That is, the three basic dimensional quantities Ci,C2 and L have been exchanged for two dimensionless parameters 
and the overall energy scale Eq. 



1 2 3 4 5 6 

<t>i 



FIG. 1: (Color online) Potential as in Eq 1221 . with its four wells. The coordinate <j>\ (target bit) runs horizontally and <j>2 
(control bit) vertically. 



III. CNOT BY ADIABATIC INVERSION 



With SQUID qubits, the logical states are the flux states of the superconducting rings ( with one qubit |0) or |I) 
with two qubits |0) |1), |1) |1) and so forth), whereas the Hamiltonian eigenstates are in general a linear combination 
of them. In case of two coupled qubits to form CNOT, each logical state of the gate will be represented by a wave 
function localized in one of the four distinct minima of the potential of Eq[22] . The four states can be labelled as 

1,2,3,4 and placed in a tableau of the kind ^) ' w ^ ere ^ e l° ca ti° ns refer to Fig. pp. That is, in the tableau the 

positions of the numbers indicate in which well of Fig 1 the state is localized while the numbers themselves indicate 
which energy eigenstate is meant. The lowest energy eigenstate is "1" , and the highest of the four first states "4" . 

A CNOT operation is defined by the conditions: A) the control bit does not change its state, and B) the target bit 
is reversed or not reversed, according to whether the control bit is |1) or |0). In the tableau representation, a physical 
embodiment of CNOT would be 



(il)-(iO « 

Condition A) on the stability of the control bit is exhibited in that no states move between the top and bottom 
row. Condition B) is realized in that the top row remains unchanged while the bottom row is "flipped" . 

Realization of operations such as Ea|2tij can be accomplished by using adiabatic processes, thanks to the "no level 
crossing" behavior of adiabatic evolution. The no-crossing property assures that a state initially in the first, or second, 
or third,.... energy level will end up in the first, or second, or third,... energy level after the adiabatic evolution, while 
at the same time the logical state associated with the level may be changing. 

One can proceed as follows: we first search for an initial Hamiltonian whose variable external parameters (</>f , 0| ) 
are adjusted to give the first four energy eigenstates localized in the four different minima of fig. 1. Then, we search 
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112=0 . 0005 




-0.01 -0.005 0.005 0.01 

^ext 

FIG. 2: A region of the 0i xt , <J)2 Xt plane with a well denned set of wavefunctions as explained in the text. The coupling parameter 
is I12 = 0.0005, The other parameters are h = h = 1, Pi = P% = 1.19, fii = fj,% = Vo = 16.3. 



for a final Hamiltonian where another set of {4>\ xt , <f>2 Xt ), gives the tableau on the right of EqJ2H]. In this procedure we 
need only to study the stationary Schroedinger equation at first. After having determined some suitable parameter 
sets, we shall also study the full time-dependent Schroedinger equation to determine what sweep speed is "slow" in 
order to guarantee adiabatic behavior. 

For the present paper we defer a detailed discussion of the phases the states acquire during the time evolution. These 
phases contain dynamic and geometric contributions pH and they themselves can be used for quantum information 
processes [22] . For a given set of parameters and sweep conditions the phases can of course be calculated explicitly 
and in no way affect the general applicability of the results and procedures presented here. Calculations of phases 
will be reported upon in future work. 

We obtain!^ the switching behavior according to Eqj2H| by fixing the control bias (f>2 Xt at a relatively high value, 
while an adiabatic sweep of the target bias 0f x * occurs. This is a generalization of a NOT 16] on 0i. The presence of 
the ?i2 coupling produces an extra bias on the target bit which "helps or hinders" the NOT operation. The relatively 
large bias on 02 comes from condition A): we attempt to "immobilize" the control bit despite the perturbations 
communicated by the sweep of 0f c * via £12. We therefore investigate the region \4>\ xt \ << 1 02**1- H 02 is indeed 
successfully "immobilized", it will be fixed in one of its two potential wells and can have only the values 02 ~ ±1- 
As seen by 0i, these two states amount to an extra bias which is added or subtracted to 0f x *. To linear order (since 
we take all § ext , cf>l xt small compared to 1 and 0i,02 are in the neighborhood of 1) and introducing the notation 
<kTeff — 0i ;rt ^ Hf tne potential terms involving 0i in Eq [22j become 

-2Z 1 1 0f rf -2/ 12 1 (±l) 

= -2Z 1 1 (0f t ± ^) = -21 x <h<l>Tt ff , (27) 

so that there is an effective shift in the external bias on 0i by (zb^ 2 -). According to Eq the bias condition for 
switching from one tableau to another is given by |0f xt | = l\2/h- 

By a numerical study where we find essentially exact solutions of the Schroedinger equation for our potential, we 
find that there is indeed a region of 0f xt , 0^' plane where the eigenstates of Eq |22 are well defined and behave as 
in this description. These regions are shown in the grey areas of Figs [21 and [3] . 

By well defined we mean that the expectation values of 0i, 02 are at the location of one of the wells w (±1, ±1), that 
only one of the first four levels is so localized, and finally that the dispersion of each coordinate y/ < <fi 2 > — < > 2 
is small compared to this expectation value of 0, that is small compared to 1. For the cases we present here the ratio 
of the dispersion to the expectation value was in the vicinity of 0.3. 

We find that the three different gray areas of Fig [2] and [5] have well defined tableaux as follows: 
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FIG. 3: As in Fig. 2 with coupling parameter Z12 = 0.005. 



3 4 
1 2 



4 3 
1 2 



4 3 
2 1 



(28) 



for the intermediate gray region (left), the darkest region (center) and the light gray region (right), respectively. On 
the order of 10-20 points were used to determine the boundaries in the figures. These tableaux fit with the description 
arrived at in the "immobilization" model, where either the top or bottom row inverts as we go from the central region 
to large |^>f |. Hence a sweep from the central region to the right region will produce the desired mapping of Eq |26j . 
Similarly sweeping from the right region to the center and from the left region to the central region and vice- versa can 
also serve as realizations, differing simply in the assignment of (0,1) for the bits or the names (1,2) for the SQUIDs. 
Note that the switch between tableaux occurs quite close to \<f>i xt \ = \li2/h\, as predicted by Eq |27| . 

Concerning the time dependent problem, the important time scale in the present context is T ac ii a b, the shortest 
time in which an operation can be performed adiabatically. This time is relevant with respect to decoherence and 
relaxation effects, since the operation must take place in a time short compared to decoherence and relaxation times. 
An estimate for T a diab for the NOT operation [lfj gives 



Tadiab — £U 



2 

tunnel 



(29) 



where e is the asymmetry of the potential and 

^tunnel — ^~osc the inverse tunneling energy or oscillation time between 
the two states. Since here we also perform a kind of NOT, we expect a similar relation to hold, where Wtunnei or t~^ c 
is the smallest level splitting during the adiabatic passage and e may be read off as the energy shift of the wells at the 
beginning and end of our sweep. A set of reasonable SQUID parameters for CNOT are L\ = 300pH, L2 = 28QpH, 



L\2 = l.8pH, C\ = C2 = O.lpF and If = /| = 1-45/xA Since in frequency units E 



1000 GHZ these 



y/L/ P H C/pF 

values give Eq w 185 GHZ and r a diab ~ 2.7 • 10 _9 s. Preliminary results are in agreement with this estimate [Ji 
and the question will be studied in more detail in further numerical work. 



IV. IMPLEMENTATION 



A. Dissipation 



As with all discussions of quantum computation, the important open questions concern dissipative effects. These 
include the decoherence time Td ec , or its inverse the decoherence rate D, as well as relaxation processes , represented 
by a time r re i ax . 
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FIG. 4: The probability of inversion as a function of sweep time. The inversion probability is one until the sweep time 
approaches Tdec, then it falls to about 0.5. When T aweep is increased toward T re iax and greater, the system always ends in the 
lowest state, so the inversion probability is one if the initial state is the ground state, while it falls to zero for starting in the 
excited state. SQUID parameters used for the simulation were /3l=1.2, the loop inductance L=400pH, the junction capacitance 
C=0.1pF, the effective resistance R=lMfi and the temperature T=40mK. 



These will affect adiabatic processes as we can see by examining the inversion or NOT process for a single qubit. 
With an increasing loss of phase coherence as caused by D, we expect the situation to become more and more 
"classical" and finally, when D is very large, the inversion is inhibited |Tt| . Furthermore, after the inversion is 
completed there will be some tendency of the upper energy state to relax to the lower energy state. 

We illustrate these effects in Fig. 0] where we show the probablity for an inversion as a function of sweep time 
Tsweeep- Note that the times involved are much longer than the above estimate for the adiabatic time of 2.7 • 10 _9 s, 
so the adiabatic condition should be satisfied. Starting at the left, we see that the probability of inversion is one until 
Tsweepi approaches Tdec- For r sweep much longer than Tdec it falls to 0.5, the "decohered" limit. For longer times, as 
Tsweep approaches and passes T re i ax , the final result depends on which state we started with. Evidently in this limit 
we will always end up in the ground state, so the inversion probability is one if we started with the ground state (solid 
line) and zero if we started in the excited state (dotted line). 

The times used were Tdec = 0.2//S and T re i ax = 0.2ms. The early part of the curve was found using the plots in Fig 
5 of the second of ref^(| with this value of Td ec , while the relaxation effects were found from the calculations shown 
in Fig 7 of this reference. The value Td ec used was based on the estimate D = T/Re 2 ^(|, with R=1MS1 at T=40mK. 

It is evident, as exemplified by the latter formula, that a necessary and-less evident- perhaps also a sufficient 
condition for quantum behavior of the system is a low classical dissipation. In the present context this implies a 
large value of R. Therefore we present some experimental data on this point, collected in the thermal regime for 
superconducting devices based on SQUIDs (24[. 

In order to evaluate the dissipation of our system, we measured the transitions between adjacent flux states of the 
rf SQUID as a function of the external flux cf> ext . In the absence of noise, the escape from the metastable well would 
occur at a critical value of the external flux <f> c . Thermal noise induces transitions at random values of <jf xt smaller 
than <f) c , whose probability distribution, namely V{4> ext ) 1 was measured by standard "time fly technique", as explained 
in the third reference listed in 

By fitting the data for P vs <p ext with Kramers theory [23 in the extremely-low-damping limit, with L,I C and C 
independently measured, we can obtain the effective resistance R. We introduce a dimensionless parameter Q ~ R, 
defined as Q=woRC, where loq is the small oscillation frequency. With decreasing temperature Q is large, showing 
a small dissipation at low temperatures. The plot of Fig. jS] indicates an exponential increase of R, and using the 
measured parameters we obtain R=22kf2 at T=2.9K. This is encouraging if the D — T/Re 2 estimate is correct, where 
we needed !Mf2 at T=40mK. The exponential fit shows that the effective resistance R is determined by tunneling of 
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B=(15.45±0.85)K 
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FIG. 5: Q factor as a function of the inverse of the temperature obtained by switching flux measurements |24|. Q increases 
exponentially with decreasing temperature, following the exponential law Q = Qoe B ^ T , with Qo = 35.37 and B=15.45K. This 
behaviour shows a strongly decreasing dissipation with temperature, and that the dissipation mechanism is essentially due to 
the tunneling of thermally activated quasiparticles in the Josephson junction. 



thermally activated quasiparticles, as expected when the external noise has been filtered out and only the intrinsic 
dissipation acts. 



B. System Design 

On-chip integrated dc-SQUIDs can be used to read out the flux states while the manipulation of the super- 
conducting qubits can be controlled by on-chip superconducting electronics. The coupling between the probe and the 
read-out system could be through a superconducting transformer. While the intrinsic dissipation can be considerably 
reduced at low temperatures as explained above 0], |lf)j | a major cause of difficulty can be the spurious interaction of 
the qubit with readout and control devices. 

In a good design, the readout device could be turned off during the manipulation. The technology to control 
this coupling may need to be developed, but could be reasonably provided by stacked junctions |20l |. or small double- 
junctions loop Jl9|, interrupting the coupling transformers. The coupling is then controlled by external current signals. 
A further possibility is the development of a fast switch using simple single flux quantum circuitry for switching the 
interaction on and off |26[ . 

In conclusion we have presented a general method for studying the adiabatic evolution of a Hamiltonian describing 
a multi-qubit system, controlled by varying external parameters. Detailed calculations were provided for a two-qubit 
Hamiltonian, whose eigenstates can be used as logical states for a quantum CNOT gate. From the numerical analysis 
of the stationary Schroedinger equation we obtained sets of parameters suitable to perform a CNOT operation, and 
indicated how a time-dependent study determines the limits for adiabatic evolution. Specializing to a definite physical 
system involving SQUIDs, we identified resaonable values of the parameters, estimated effects due to dissipation and 
considered some points of system design. 
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